General Theory of the Quantum Kicked Rotator. I 



Tao Ma 

Department of Modern Physics, University of Science and Technology of China, Hefei, 

(Dated: February 2, 2008) 

This is the first of a series of two papers. We discuss some basic problems of the quantum kicked 
rotator (QKR) and review some important results in the literature. We point out the flaws in the 
inverse Cayley transform method to prove dynamic localization. When t/2tt, where r is the kick 
period, is very close to a rational number, the localization length is larger than the typical localization 
length. We analytically prove anomalous localization and confirm it by numerical calculations. We 
point out open problems that need further work. 
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PHYSICAL EXPLANATIONS OF QUANTUM 
KICKED ROTATOR 



Nearly thirty years ago, QKR is first studied by G. 
Casati, B. V. Chirikov, F. M. Izraelev and J. Ford 
They discovered by numerical calculations when the kick 
period r is the product of an irrational number and 2tt, 
the rotator localizes in the momentum space. Later, S. 
Fishman, D. R. Grcmpcl and R. E. Prange explained 
the localization by transforming QKR into an Anderson 
localization problem 2|. It does not seem necessary to 
discuss the basic problems of QKR again. But many 
problems are impossible to be solved by present methods. 
The paper is both a review and a problem list. 

The Hamiltonian of QKR is defined as [1| 



H = H + V(9) ]T 8(t - tit) 

n=l 
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which describes a particle restricted to a ring with free 
Hamiltonian Hq = — \^ 2 -§q2, and periodically kicked by 
a homogeneous electric field parallel with the ring plane. 
The parameters r is the kick period and k is the kick 
strength. In QKR, the interaction between the electric 
field and the rotator is V{6) = —kcosd. But in this 
paper sometimes we discuss a more general form of the 
interaction V(9). In the basis {\m) = -±=e~ ime }, the 
Hamiltonian can be written as 
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The following kick system has the same Floquet oper- 



ator as the Hamiltonian in Eq. (2) 
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The physical explanation of the Hamiltonian in Eq. 
(3) is a particle on a one dimensional lattice. \m) is 
explained as the m— th Wannier state or m— th site. 

The term - E™=-oo(|l™><™ + W m + 1>H) * 
the hopping matrix or the kinetic energy. The parti- 
cle is periodically delta-kicked by a harmonic potential 

En=i E™=-oo T^\m)(m\S(t - nk). Now k is the kick 

period or the free diffusion time of the particle, and the 

2 

potential is defined as 

Eq. (3) is discovered when we solve the quantum 
kicked linear rotator (QKLR) 0, 0| ■ QKLR always lo- 
calizes except when r = 27r x Integer, where Integer is 
an integer. There is Bloch oscillation phenomenon when 
a lattice is put in a homogeneous electric potential. It 
turns out QKLR localization is a general Bloch oscil- 
lation. QKLR and the Bloch oscillation are the linear 
Toeplitz system in two different representations: the ro- 
tator representation and the site representation. QKR 
can also be explained in the site representation. 

The site representation of QKR or Eq. (3) can be an- 
other implementation in the laboratory compared with 
the usual implementation Eq. (2) 0, @. The classical 
correspondence of the site explanation is a periodically 
delta kicked classical random walker on a lattice by the 
harmonic potential (CKRW). In the rotator representa- 
tion, the classical correspondence is the standard map. 
We think CKRW is also very interesting just as the stan- 
dard map. For example, does CKRW localize when the 
kick period r is the product of an irrational number and 
2-7T? Gong et al independently gave an equation similar 
to Eq. (3) and proposed QKR can be implemented as 
a Heisenberg spin chain subjected to a parabolic kicking 
magnetic field @, [f| • See Q for detailed information. 
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The general kick system can be defined as 

oo 

H = H + V^5{t-nT). (4) 



The Floquet operator is F = e 



tH t 



We can al- 



ways treat V as the unperturbed Hamiltonian and Hq the 
perturbation just we reexplain Eq. (2) as Eq. (3). The 
essence of delta kicked system including QKR is two dif- 
ferent Hamiltonians acting on the Hilbert space in turn. 

The following Hamiltonian has the same Floquet op- 
erator of Eq. (1). 



H{t) = 



-¥ 



2 d 2 
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< tMod(r 
t < tMod(r - 



V k) < t 
k) <t- 



(5) 



where tMod(r + k) is the product of the fraction part of 

TTfe and T + k - 

Apparently n different Hamiltonians 

{Hq, Hi, ■ ■ ■ ,if n _i} can act on the wave func- 
tion in turn. Then the Floquet operator is just 
F = e -| ff «-i T "-i . . . e -tffiTi e — |H r _ The a b ove 

explanation also applies classically. 

In a laboratory it is impossible to implement Eq. (5) 
because Hq always exists. So the experimental imple- 
mentation is we increase the interaction strength and 
decrease the interaction time. When the interaction 
strength is very large and the interaction time is very 
small, it is almost a delta function. 

When M is very large the following system is the ex- 
perimental implementation of Eq. (1). 



H{t) 



2 n W' 



\Ho 



< tMod(r + -^) < t 
M cos 9 t < tMod(r 
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or 



H(t) 



2 n ae 



< tModir 
kM cos 9 t < tMod(r - 
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We define one period as a free evolution of period r and 
a kick, so after n periods, the QKR wave function is 



V>(nr) = r^(0), 



(11) 



where the unitary operator U (n) = F n maps the initial 
state to the state at the time nr. 



A. Meaning of dynamic localization 

Now we ask the fundamental problem of QKR. If the 
rotator is initially in the ground state, will its energy 
SreL-oo \ \ c n\ 2 n 2 increase to infinity when the time runs 
to infinity? This problem can also be described in an- 
other way, such as the wave function X)^L-oo l c ™ I ^ s nor " 
malizable. The two questions has some nuanced differ- 
ences. We do not know whether another situation could 
happen. For example, its energy is infinity, while its wave 
function is still normalizable. This situation may happen 
under some conditions, but in this paper we assume, to 
QKR, the energy finity is equivalent to the wave func- 
tion normalizability. And the situation of energy finity 
and wave function normalizability as time runs to the 
infinity is referred as dynamic localization. 

There are close analogies between QKR and some 
problems in the solid state physics. The quasienergy 
band of the Floquet operator with rational Z- is anal- 
ogous to the energy band of the Bloch operator. The 
irrational case Floquet operator has been transformed 
into an Anderson localization problem in (2j. But there 
is one fundamental difference. An Bloch electron can be 
in the eigenstate of the Bloch operator, while QKR can 
never be in a Floquet eigenstate, because a Bloch eigen- 
state has a finite energy, while in the case of QKR the 
a Floquet eigenstate has an infinite energy. So the rota- 
tor state will always be a superposition of many Floquet 
eigenstates. The absolute values of the coefficients of the 
eigenstates in the superposition never change, and what 
changes is only the relative phases between different Flo- 
quet eigenstates. 



II. FLAWS OF THE INVERSE CAYLEY 
TRANSFORM METHOD 

We study the Schrodinger Equation 

ih^i(>(6,t) = Hip(0,t). (8) 

The dynamics of a free evolution of period r and a kick 
is described by the Floquet operator F. 

F = e-i y W e -*^ (9) 
The matrix elements of F are given by 

F nm = (n\F\m) = exp(-^r^-)i m -" J w _ m (|). (10) 



B. Inverse Cay ley transform and dynamic 
localization 

One way to study QKR is to find all the eigenvalues 
and eigenstates of the Floquet operator. If the eigen- 
states are extended, dynamic derealization happens and 
vice versa. Milek et al has formally proved absolutely 
continuous spectra of the Floquet operator imply der- 
ealization. When the Floquet operator has singular con- 
tinuous spectra, the usual dynamic derealization does 
not happen @ . Can the Floquet operator have extended 
eigenstates and singular continuous spectra under some 
conditions, for example, when ^ is a Liouville number? 
Jitomirskaya et al proved almost Mathieu equation has 
singular continuous spectra under some conditions p^| . 
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Yet we assume localized eigenstate is equivalent to the 
dynamic localization to QKR. 

One way to find all the eigenstates of the Floquet oper- 
ator is to transform the Floquet eigenstate problem into 
time independent Hermitian operator eigenstate prob- 
lem. Especially when ^ is irrational it is transformed 
into an Anderson problem (2j. The Floquet eigenstate 
problem is 



,-iA, 



(12) 



A Hermitian operator U can be transformed into a 
unitary operator O by Cayley transform O = A 
unitary operator O can be transformed into a Hermi- 
tian operator by inverse Cayley transform U — ijr%- 
In complex analysis, Cayley transform, which is a linear 
fractional transform, maps the real line to the unit circle 
on the complex plane and the inverse Cayley transform 
maps the unit circle to the real line. Cayley transform 
and its inverse transform of both the operator form and 
the complex number form has many applications in dif- 
ferent fields of mathematics. 

Substitute 



1 + iU(0) 

i - iu{oy 



(13) 



where U{6) = -tan(^) 
(12). We get 



tan(-|cos#), into Eq. 



T m u\. 



■5> 



= 0, 



(14) 



where u\ = (1 + e'^ A H °')4>x, U r — U- r is the Fourier 
coefficient of U{6), and T m = tan x ~ T 2 m ~ ■ QKR is trans- 
formed into a particle moving on a periodic or non- 
periodic lattice depending on r. The evolution of QKR 
has no direct relation with the new particle on a lat- 
tice. But the eigenstates of the Floquet operator and the 
eigenstates of the new formed Hermitian operator are re- 
lated by Ux = (1 + e l ( x ~ H °))(t)\. If 4>\ is extended or 
localized, then u\ is also extended or localized and vice 
versa. 

Wc define A as 



T_ 2 Ut U 2 U 3 

u-i r_i Ux u 2 

U- 2 U-t T Q Ui 

Us U-2 U-i Ti 

■■■ U-3 U-2 U- 



Ui 
T 2 



\ 



(15) 



Then 



Au x + U u x = 



(16) 



One basic observation is the diagonal matrix elements 
T m of A are pseudorandom numbers if ^ is irrational. 
T m is a Cauchy distribution. Because U(6) looks quite 



regular here, now we invoke Anderson's result [llj. u\ is 
localized. The above analysis gives an analytical proof of 
dynamic localization discovered in [l|. Fishman, Grem- 
pel, and Prange's inspiring result connects the fields of 
quantum chaos and Anderson localization in disordered 
matters. It points out the physical mechanism of dy- 
namic localization of QKR. Inverse Cayley transform 
seems to be the only method to transform the Floquet 
operator eigenstate problem to the disordered Hamilto- 
nian eigenstate problem. 



C. Flaws in the inverse Cayley transform method 

Now we check the above argument step by step. 
The equidistribution of the sequence {an 2 Modl} n when 
a is irrational can be proved from the equidistribu- 
tion of sequence {an} n and van der Corput's Theo- 
rem. While the equidistribution of sequence {anModl} n 
can be proved by Weyl criterion. So all the sequences 
{an m Modl} n , where is to is a positive integer and a 
irrational, are equidistributed between and 1. Other 
equidistribution examples include {x n Modl} n , where 
x > 1. The sequence of all multiples of a by all 
prime numbers {2a, 3a, 5a, 7a, 11a, . . . Modi} is stud- 
ied in [12j as the quantum kicked prime number ro- 
tator (QKPR). It localizes when a is irrational be- 
cause {2a, 3a, 5a, 7a, 11a, . . . Modi} are equidistributed 
between and 1. QKPR also localizes when a is rational 
such as | because now it is analogous to a generalized 
kicked dimer model. 

The number theoretic property of a number sequence 
is used in Shor's algorithm to find prime factors of a 
composite number S. If S is a large number, and we 
define a rotator with energy levels {p n ModS} n , where 
n runs from 1 to S and 1 < p < S. If S is a prime 
number, dynamic localization (To a finite Hilbert space, 
this is not very rigorous) also happens. If S is a composite 
number, dynamic localization generally does not happen. 
Although we can not find the exact prime factors of S, 
at least we can judge whether S is prime. 

There is one inherent paradox (Local pseudorandom- 
ness paradox) of localization caused by number theoretic 
randomness. On the one hand, the randomness of the 
sequence {^n 2 } n with irrational j-, is only meaning- 
ful when n goes to infinity. On the other hand, if the 
wave is localized, how could the wave feel the random- 
ness of j^n 2 of very large n, where there is randomness? 
This can be said in another language, how can a remote 
(global) randomness influence a localized wave or how 
can a localized wave feel a remote (global) randomness? 
In Anderson localization, every random variable is inde- 
pendent. So in the place where the wave is localized, the 
randomness is real randomness. But in number theoretic 
randomness, there are strong correlations between the 
local pseudorandomness. For example, the QKR wave 
function is localized at several localization lengths. Is the 
pseudorandomness at several localization lengths randon- 
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mness or even pseudorandomness? Let's take k — 1. So 
the wave function is localized between several basis vec- 
tors. For example, they are from | — 10) to 1 10) . Is the 
sequence {-^n 2 } n from n = —10 to 10 pseudorandom- 
ness? 

Randomness can be classified into two categories, the 
statistical randomness and pseudorandomness caused by 
chaos or number theoretic origin. The local pseudoran- 
domness paradox is not a problem which only concerns 
QKR or Eq. (16). It is also not a new question. L. Boltz- 
mann tried to found statistical mechanics on the crgodic 
theory. A recent example is Bohigas-Giannoni-Schmit 
conjecture (l3| . 

The first problem of the inverse Cayley transform 
method to prove localization is to some irrational , the 
sequence {^n 2 } n may be not random enough to cause 
localization. This is the pseudorandomness problem that 
has been discussed by many authors, such as P.[l3.[l5j. 
In [T3 |. we give an example of irrational 4= which is not 
irrational enough to cause localization. 

Second, Eq. (16) is not really an Anderson localization 
problem. In Anderson localization, we know the Hamilto- 
nian, and we do not know the eigenvalues and eigenstates. 
While in Eq. (16), we know the eigenvalue beforehand, 
which is just —Uq. While the unknown A is contained in 
the diagonal matrix elements T m . This is not an Ander- 
son localization problem. It is rather an inverse Ander- 
son localization problem. We know one special eigenvalue 
and we need to find what diagonal "disorders" satisfies 
Eq. (16). Surely there are infinitely many diagonal "dis- 
orders" satisfies Eq. (16). Eq. (16) has the form of an 
eigenvalue problem, while in fact it is a nonlinear equa- 
tion of A or e~ lX . One may say to a general irrational 
whatever A is, the sequence T m is random enough to 
guarantee the eigenstate localized. In Anderson's formu- 
lation of localization, the diagonal matrix elements are 
required to be independent random variables (ll| . No 
independent random variables can satisfy an equation, in 
which these independent random variables are dependent 
variables. We use the terms independent and dependent 
variables in their original sense of Statistics. 

Generally an infinite matrix has infinite eigenvalues 
and eigenstates. For example, we find a A which satisfies 
Eq. (17). Then the operator A defined in (16) surely has 
many other eigenvalues besides — Uq. Are other eigen- 
states of the infinite matrix with definite A also localized? 

The third problem is U{9) — — tan(| cos 9) is regular 
only if k < it. When k > tt, U(9) is discontinuous at some 
points. The first discontinuous point is 9 = arccos 
At this point U(9) is oo. In [ll|, the interaction is re- 
quired to be "falling off as the distance r — > oo faster 
than 1/r 3 ". In [ll| . r is the distance between two lattice 
sites, while r in U r is the index. But both two express 
the interaction strength between different sites. Does U r 
fall faster than 1/r 3 ? This question can be answered 
in two ways. First, if U r < where c is a constant, 
then U{9) = J2r U r e lr9 is convergent at every point of 
9. So U r can not fall faster than 1/r 3 . Second, when 



k > 7r, U(9) is not square integrable. We only need to 
check whether U{9) 2 is integrable in the small domain 
arccos f — 6 < 6 < arccos f + S, where S is a small num- 
ber. If k > 7r, in the domain, %V(6) can be approximated 
by | +c9, where c is the derivative of | V(9) at the point 
9 = arccos | . But tan( : | + c9) is not square integrable 
in the domain. If k = tt, tan(ifcV (9)) 2 is approximated 
by (tan(f(l - f ))) 2 = (cot(fx 2 )) 2 « l/(sin^) 2 » 
1/(t~) 2 - It is not square integrable either. So when 
k > tt, U(9) is not square integrable. J2 r U 2 is divergent 
due to Parseval's identity. The worst estimate is U r falls 
slower than \j\fr. The Fourier expansion U r of U(9) 
does not exist at all, when k = tt x Odd, where Odd is 
an odd positive integer such as 1,3,5,.... 

So when k > n, we can not invoke Anderson's result 
to prove the vectors u\ are localized in Eq. (16). This 
is a fatal flaw of the inverse Cayley transform method to 
prove dynamic localization. In one dimension it is easy 
to localize [17]. So the slowly falling U r may not destroy 
localization. But in the conventional theory of condensed 
matters, the strongest long range force is Coulomb force 
oc ~. QKR with a general irrational number ^ seems 
to always localize however large k is. We think the real 
meaning of the breakdown of the inverse Cayley trans- 
form method is there may be a localization-delocalization 
transition when k increases from k < tt to k > tt and 
k = 7T is the critical point to many kicked systems. In 
the QKPR of [12], there is apparently a localization- 
delocalization transition in the rational case r = 
When k = 1 the prime number rotator localizes, while 
when k = 5 it delocalizcs. 

There are several flaws of the inverse Cayley trans- 
form method. Why not develop an independent local- 
ization theory concerning Floquet operators or unitary 
operators? And the new theory which treats time de- 
pendent problems (unitary operators) is parallel with 
the Anderson localization theory which treats time in- 
dependent problems (Hermitian operators). This does 
not seem to be an easy task. The first difficulty is there 
is not a unitary perturbation theory concerning unitary 
operators, while to Hermitian operators there are many 
perturbation methods. We will return to the problem of 
developing a dynamic localization theory of the Floquet 
operator in the second paper. 



III. ANOMALOUS LOCALIZATION 

The inverse Cayley transform method provides a phys- 
ical picture, because of diagonal pseudorandomness, the 
eigenstate tends to be localized. But there is another 
path of the development of QKR theory, which gives a 
different picture. Casati and Guarneri proposed there is 
a non-empty set of irrational , to which dynamic local- 
ization can not happen [3]. We must resolve these two 
conflicting views. Fishman et al also pointed out to the 
Liouville number things are very delicate and they 
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excluded the Liouville number ^ in their theory [2|. 

We have to discriminate two kinds of different "der- 
ealizations" . To general irrational jr , the localization 
length l T is estimated to be aD (ill [l9j . where D is 

,2 

the classical diffusion constant And the factor a is 
estimated to be \ in [2(| Hl| . The first kind of "dereal- 
ization" is aD < l T < oo, which is actually localization 
with an anomalous localization length. This is referred 
as anomalous localization. The second derealization is 
If oo . The anomalous localization can be seen from ex- 
amples such as -£- = i + n^ni or the continued fraction 

1 47T O ZX1U 

j- = (0; 3, 10 m , 1, 1, 1, • • •), where m is a large number in 
both cases. 

But we do not know all the irrational numbers ^~ 
with anomalous localization or derealization. We do 
not know to a special J- whether anomalous localiza- 
tion or derealization depends on k. For example, when 
k < 7r, derealization does not happen to a t, while when 
k > 7r it does. We do not know whether both cases 
can happen to general Liouville numbers, such as Liou- 
ville constant L — J^JLi 10~-' ! . The continued fractional 
of L is {0; 9, 11, 99, 1, 10, 9, 999999999999, . . .}. The n-th 
incrementally largest term consisting only of 9s occurs 
precisely at position 2 n — 1, and this term consists of 
(n- l)n! 9s [ll- 

Although we have constructed an irrational ^ with 
derealization in [l6| . There is not one specific irrational 
^ that has been numerically calculated to be delocal- 
ized. It is difficult to do a numerical simulation to these 
special irrational numbers. It is the periodic structure of 
the Floquet operator which causes derealization. So a 
successful calculation must preserve the periodic struc- 
ture. To the rational 4- — - case, the truncated Hilbcrt 

Z7T q 

space is at least as large as several qs. This makes it 
very difficult to calculate the rational cases with large 
denominators. 



A. A quantum analogy of Lyapunov exponent 
equation 

There are two different Floquet operators F and F 1 . 
We calculate the difference between U(N) — F N and 
U(N)' = F' N . We define SU(N) = U(N)' - U(N), SF = 
F'-F and 5t = t' - t. 

SU{N) = F'U{N - 1)' - FU(N - 1) 
= F'U{N - 1)' - F'U(N - 1) 
+ F'U(N - 1) - FU(N - 1) 
= F'(U(N - 1)' - U(N - 1)) + SFU(N - 1) ^ 
= F'8U{N - 1) + SFU(N - 1) 

N 

= J2u(N-j)'6FU(j-l). 
i=i 

In the last step, the mathematical induction method is 
used. Eq. (17) describes the divergence between unitary 



operators because of small difference SF. This is a quan- 
tum analogy of the Lyapunov exponent equation. The 
classical Lyapunov exponent equation describes the di- 
vergence between two orbits in the phase space because 
of small difference between the initial conditions. The 
stability of quantum systems is indicated by their sensi- 
tivity to the perturbation of the Hamiltonian [23| . We 
think the sensitivity of the Floquet operator to a per- 
turbation is another criteria to judge whether chaos is 
relevant to a quantum system. In other words, the sta- 
bility of a periodically driven quantum system is decided 
by the sensitivity of the Floquet operator to perturba- 
tion. There are two kinds of sensitive perturbation to F. 
The first kind is the phase sensitivity. The perturbation 
rotates the phase of some matrix elements by an angle 
7T. The second kind is the perturbation greatly modifies 
the module of some matrix elements. In QKR, SF is not 
sensitive to Sk or 5h. It is sensitive to St. We consider 
St = t' — t. The matrix element 

JV 

5C7(iV) nm = ^^(C/(JV-i)'e-^) n! 

Hi ( l8 ) 

xie-^'-e-V^UU-l)^. 

For simplicity, we now assume n — 0. The summation of 
I of U(j — l)i m is effective in a finite bandwidth, which is 
at most from -jk to jk. | e -i/2^V_ e -i/2 4 ; 2 r| ^ \/21 2 5t. 
The row (U(N - j)'e~ iV ) n i and the column U(j - l) im 
are unit vectors. From Cauchy-Schwarz inequality, 

| J2(U(N - 3)'e- W ) n i{e- 1/2il2T ' - e~^)U{j - l) im | 2 
i 

<CE\(U(N-3)>e-* v U 2 ) 
i 

x{ J2l(e-^'-e-^)U(j-l) lm \ 2 ) 
i 

= ^ |(e -l/2^r'_ e -l/2«V )[/(i _ lW|2 
I 

<]T|l/2Z 2 <M/(j-l) im | 2 
( 

= ]ri/4Z 4 <ST 2 |[/(j-l) im | 2 
l 

< l/AijkfSi 2 , 

(19) 

where the equal sign = in the last step comes from when 
\U(j — 1); TO | 2 concentrates on the boundary jk of the 
bandwidth, which is actually impossible. 

| Y,(U(N - j) / e- <V k(e- 1 / aaV - 6 " 1/m2t ) 

l (20) 

U(j - l) lm \ < l/2(jk) 2 ST. 
After the summation of j in Eq. (18), 

\SU(N) nm \ < l/6N 3 k 2 ST. (21) 
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Before the time (til- 
first got by Casati et 



1/3 



JSU(N) nm \ < e. Eq. 
(14| . In fact, generally 



\SU(N) nn \ < l/3iW<J- 



(21) is 



(22) 



To the general rational number the diffusion speed 
is far slower than the case rModAn = 0. So the effective 
bandwidth of U(j — 1) is far smaller than jk. The diffu- 
sion speed v of the case r = is k. To a rational the 
diffusion speed is v T . Then the bandwidth of U(j — 1) is 
v T j. So 



|<ft/(iV)„ ro | < 1/6N^St 



(23) 



When t/r < (-^p) 1 ' 3 , \SU(N) nm \ < e. We can define 
(^jfp) 1 ^ 3 T as the divergence time, after which the U(N) 
and U(N)' are significantly different from each other. v T 
is quite small except some strongly resonant cases. To the 
resonant case = - with k = 1, we estimate v T can be 
as small as cie _C2<3 , where C\ and C2 are two constants, 
the relation between which and q and k is unclear at 
present [l6j . 

In the derivation of Eq. (21), the summation of I 
is from — jk to jk. This is the upper limit of diffu- 
sion with the fastest diffusion speed, which happens only 
when r = or 47r. The more exact value of the band- 
width is %jk from the property of Bessel function, where 

e = 2.71828 • • • . So a factor ^ = 1-84726 • • • should be 
added to Eq. (21) and (23). We do not consider the 
factor, which will cause qualitative differences to our dis- 
cussion. 

The above calculation crucially depends on the almost 
band structure of F, F', U(N) and U{N)' . A band ma- 
trix can not be a unitary matrix. The different columns 
or rows of a unitary matrix are orthonormal with each 
other. A band matrix can never has its columns or rows 
orthonormal with each other. For example, let b denote 
the end of one row of a band matrix. Then another row 
of the band matrix has its start at b. These these two 
rows can not be orthonormal with each other. Another 
thing worth notice is e~~ lV is a Toeplitz matrix. U in Eq. 
(14) is also Toeplitz. Toeplitz matrices are quite rigid. 
The quantum kicked harmonic oscillator [24| does not 
have such a nice property. Casati et al tried to explain 
dynamic localization from the perspective of band ran- 
dom matrix [25j . The Floquet operator of QKR is neither 
band nor random. When k > n, from the perspective of 
Anderson localization, u\ should be delocalized, but it is 
still localized. The mechanism of the localization of Eq. 
(16) and QKR is even stronger than the mechanism of 
Anderson localization. It is this mechanism that QKR 
theorists need to find. 

We can get a more exact estimate than Eq. (21) under 
some conditions. In Eq. (19), we use Cauchy-Schwarz 
inequality to drop U{N — j)' . We can also drop U(j—1). 
A better estimation is to drop the unitary matrix with 
a larger bandwidth. For example, if -Z- and ^ are all 
irrational numbers, the summation of I is restricted to 



the localization length l T . 

\SU(N) nm \ < 1/6NI*6t. 



(24) 



So the divergence time between different irrational are 
far larger than the divergence time between rational and 
irrational ■S-. 

Eq. (24) also holds when a rational Z- is close to a 

general irrational number. For example \-%= — ^f 1 1 < e, 
where is rational and e small. 

2.1T 



Analytical and numerical proof of anomalous 
localization 



Now we estimate the anomalous localization length. 
From Eq. (23), before the divergence time 



VtOT 



(25) 



QKRs of r and r' will not diverge from each other much. 
We assume ^ is rational and ^ irrational. So QKR of 
t' will at least diffuse to the length (^f^) 1 / 3 . We assume 
e = 1. Then the localization length l T i is at least (^) 1/,s . 
Another method is from Eq. (24), if the divergence time 
is Nt, then N = jrj^:- And l T i — Nv T . We get the same 

result. So 



lr> = (^) 1/3 - 
or 

We assume v T (k) = v T (k = \)k. Then 



lr- = ( 



6v T (k = l)fc 1/3 



St 



(26) 



(27) 



The anomalous localization is significant when v T is 
large, and St small. v T is large when r is close to strong 
resonances and k large. In our numerical calculation, we 
take is the sum of ^ and a small number, r = 4tt^ has 

4-n- 3 3 

the fastest delocalization speed except r = or An. In 
our calculations, K = 1, k = 1 and the truncated Hilbcrt 
space is from | — 500) to |500). 

In FIG. 1 and 2, r = 4?r x [0; 3, 100, 1, 1, 1, ■ • ■], 
where the continued fraction [0; 3, 100, 1, 1, 1, ■ ■ •] w 1/3— 
0.00110064. The localization length is around 6 from 
FIG. 1 and 2. 

In FIG. 3 and 4, r = 4tt x [0; 3, 1000,1, 1,1, • • •], where 
[0; 3, 1000, 1,1,1, •••] w 1/3 - 0.00011101. The localiza- 
tion length is around 15 from FIG. 3 and 4. Notice there 
is staircase like structure in FIG. 3. It is very rough. We 
do not know what they are. 

In FIG. 5 and 6, r = 4vrx [0; 3, 10000, 1, 1, 1, ■ • •], where 
[0;3, 10000, 1,1,1, •■•] w 1/3-0.00001111. The localiza- 
tion length is around 36 from FIG. 5 and 6. Note far 
outside the localization length, the wave function is not 
smooth. 

It is difficult to tell the exact localization length, which 
seems to change at different times. Compare FIG. 2, 4, 
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FIG. 1: QKR wave function at different time. N indicates at time 2 N r. For example, TV = 10 means at time 2 10 r. 
n is the n-th basis \n) = —^=e ine and c n is the base- 10 logarithm of the absolute value of the wave function on the 
n-th basis, k = 1, r = 4ir X [0; 3, 100, 1, 1, 1, ■ • ■ ]. The initial state is |0). 
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FIG. 2: QKR wave function at different time. FIG. 1 and FIG. 2 are the same calculation. 
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FIG. 3: QKR wave function at different time, fe = 1, r = 4n X [0; 3, 1000, 1, 1, 1, ■ ■ ■ ]. 



The initial state is |0). 
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FIG. 4: QKR wave function at different time. FIG. 3 and FIG. 4 are the same calculation. 



and 6; the general trend of increasing localization length 
is clear. Compare these localization lengths together. 
6 3 : 15 3 : 36 3 = 1 : 15.6 : 216. The general trend l r , cx 
(fr -1 / 3 is confirmed. It seems j^tts is a monotonically 
increasing function of i5r -1 / 3 . 



C. Casati and Guarneri's argument 



In [14| , Casati and Guarncri defined a quantity 
R T (n, tp) = i Y^j=i \(Ft V'i VOI 2 to measure the recurrent 
behavior of QKR. They proposed if the Floquet opera- 
tor F T , where r is the kick period, has purely continuous 
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FIG. 5: QKR wave function at different time, k = 1, r = 4tt X [0; 3, 10000, 1, 1, 1, • • • ]. The initial state is |0). 
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FIG. 6: QKR wave function at different time. FIG. 5 and FIG. 6 are the same calculation. 



spectrum, then when n — > oo, R T (n,ip) — > 0. It follows 
from Eq. (21) that \R r >{n,i/>) - R T (n,^)\ < 2 7 n 3 . Then 
starts Casati and Guarneri's argument. 

Rr>{N n ,i>) = \R T ,{N n ,TP) - R Tn {N n ,^)\ + R Tn {N ni iP) 
= 7k' - r„\N n + Rr n (N n ,ip) 

(28) 

Find the sequences r„ and N n . They satisfy the condi- 
tions, when n — > oo, \t' — T n \N n — * (Condition 1) and 
R Tn (N n ,ip) — > (Condition 2) at the same time. Then 
Rr'{N n ,ip) — ► 0. In this way they argue there are some 
irrational ^s, which have a purely continuous measure. 

For a long time, Casati and Guarneri's surprising ar- 
gument puzzles QKR theorists [l5|. One gap in the ar- 
gument is whether when n — > oo, Condition 1 and 2 are 
true at the same time. r„ is not a static rational num- 
ber, so R Tn (N n , ip) — > is true only if N n is large. (See 
[l6| for why N n needs to be large.) But now we can not 
guarantee Condition 1. Is there at least one irrational 
number that satisfies Condition 1 and 2? The answer 
to the question is not a priori true. 

In [lfj], we construct one irrational with delocaliza- 
tion. Our construction depends on one assumption and 
three facts. We assume the rational case has a nearly 
constant diffusion speed v T from t = to t — > oo. At 
least the rational case diffuses with a nearly constant 
speed v T after a threshold time tLD- This is the mo- 
mentum linear diffusion assumption, which is different 
from classical diffusion. First, there is an upper limit 
of momentum diffusion speed. Second, in all 77- = - 

with different p, the case = - has the lowest diffusion 

1 ' Z1T q 

speed. Third, Eq. (21) connects rational and irrational 



cases. The constructed is not a general Liouville num- 
ber. To construct more irrational we have to exactly 
know the diffusion speeds of QKR with every rational ^- . 

Because some ir- = 2 are close to 2? with q' <C q, the 

27r q q' * ^ ' 

diffusion speed v T can be far faster than other x = 
where is not close to a rational number with smaller 
denominators. This can be seen from Eq. (21) before the 
divergence time ^ = 2 and ^ will not depart from each 
other very much. This is one reason why the diffusion 
speed v T and v T i is not significantly different from each 
other. It is difficult to estimate the diffusion speed of the 
general £ = f . 



D. Some problems 



For all rational numbers = | between and 2, or 



2 has the fastest diffusion speed vq = k. For X = 0, 2, 
U(N) = F N = e -iNkV(9)_ For x. = l, the QKR is 
periodic [l[ which is proved in Q. How to sort all the 
rational numbers according to their respective diffusion 
speeds? This sequence is helpful to understand QKR. Is 
this ordering dependent on kl 

For r close to strong resonances such as {|, |, |} X 4w, 
QKR will diffuse quickly before the time (p^) 173 . This 
is the reason of anomalous localization. Even if the ma- 
jority of irrational numbers (we assume the Lebesgue 
measure is 1) are localized, irrational numbers with typ- 
ical localization length can not have Lesbegue measure 
1. All the irrational numbers between, for example, 
[47r/3 — 2e,47r/3 + 2e], where e is a small number, do 
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not have the typical localization length. From Eq. (26), 



2c 



47r/3+2e 



l T i dr' = —l T i- 



4tt/3 



Air /3+2£- 



(29) 



If ^ r '=47r/3+2e is the typical localization length, the aver- 
age localization length between [47r/3,47r/3 + 2e] is ap- 
parently not the typical localization length. Is the av- 
erage localization length I = -L l T dr finite? Even 
if it is finite, it is not the usual localization length aD. 
In Eq. 29, the small domain contributes the factor |. I 
is larger than aD because of the contribution from the 

l T dr 



small domain around r = 47r/3. Does the J Q 
also contribute a factor |? I = ^aDl 

The average localization length I is very difficult to 
calculate. The question can be answered only after we 
have totally understood delocahzation of all the rational 
t^t. Surely if we did so, we would totally understand 
QKR. We guess the average localization length could be 
infinity. Is l T dr/\Ti — T2I finite, where t\ and T2 are 
boundaries of a very small domain? We guess it could be 
also infinity. I may be 00, but the Lesbegue measure of 
finite I is still almost 1. The probability distribution of 
P(l) is surely interesting. Our guess is inspired by An- 
derson. Anderson emphasized in disordered matters the 
average localization length may be infinite, nevertheless 
the Lesbegue measure of finite localization length is still 
larg e. It is the P(l) rather than I that are really relevant 
[26|. But the problem whether 7 of QKR is finite is still 
interesting. 

Although dynamic localization happens to irrational 
numbers, the key to understand localization rests on un- 
derstanding delocahzation of rational numbers, especially 
the delocahzation speeds. In the paper, we used the top- 
down method from rational cases to irrational cases to 
understand irrational cases. Now we outline the bottom- 
up method from the irrational cases to understand ratio- 
nal cases. We use Eq. (21) in reverse. For a typical irra- 
tional number with localization, 1 7? — - 1 < , where is c 

' 1 2-7T q 1 q A ' 

is a trivial constant, according to Dirichlet's approxima- 
tion theorem on diophantine approximation. So = 2 

is still localized before the time approximately g 2 / 3 . If 
we want stronger result, we have to prove there is a typ- 
ical localized irrational It? — -I < In this way 

2?r ' 1 2tt q 1 q n J 

we may increase the time to approximately q n ^ 3 . This 
is petitio principii. But the essence is both rational and 



irrational numbers are dense within each other. If typical 
irrational number with typical localization length has a 
Lebesgue measure approximately 1, why is there not one 
such number in [~ — | + ^-]? 

In the finite time of an experiment in a laboratory, 
■j— = jr + y^-, n,i , where m is large, can diffuse much faster 

47T 3 2X1U' 07 

than an rational -f- = which has a large denominator 

4-7T q ' G 

q and ^ is far from any strong resonance. Another inter- 
esting case is when r close to or 4ir. Before the diverge 
time between r and (or 4ir), QKR will diffuse with a 
large speed. This is discussed from classical perspective 
in [27] . For small time t, the energy grows quadratically. 
When t is large and not too large, the energy growth will 
stop because of big denominator of j- or 2 — ^ [3| . This 
is referred as dynamical freezing in 27j. It is difficult to 
estimate the energy linear growth [27 1 at the intermediate 
time from our discussion here. 



IV. CONCLUSIONS 

In the paper, we discuss some basic problems of QKR 
theory. We point the flaws of the proof of dynamic lo- 
calization in Fishman et aPs method 0. We emphasize 
the physical mechanism of Anderson localization can not 
totally explain dynamic localization of QKR. We empha- 
size it is necessary to understand the delocahzation of all 
the rational -£-s. In [r|, we have numerically calculated 
the delocahzation of r = jg. Yes, it delocalizes. But 
the delocahzation time is around 2 30 r or 2 45 t. In [r|, 
we constructed an irrational 77- with delocahzation. In 

27T 

the paper, We theoretically prove anomalous localization 
and numerically confirm it. These three phenomena tell 
us the QKR theory is not just localization with irrational 
^- and delocahzation with rational The whole pic- 
ture is more complete only if when we take these three 
facts into consideration. What we have touched is only 
a fraction of the QKR theory. We point out the open 
problems and hope the readers can solve them. 
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